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f^"^ ' Abstract 

I We prove that all stationary and spherical symmetric black hole solutions to theories 

■ with symmetric target spaces are integrable and we provide an explicit integration method. 
I This exact integration is based on the description of black hole solutions as geodesic curves 
i on the moduli space of the theory when reduced over the time-like direction. These 

. ^ I geodesic equations of motion can be rewritten as a specific Lax pair equation for which 

■ mathematicians have provided the integration algorithms when the initial conditions are 
^ I described by a diagonalizable Lax matrix. On the other hand, solutions described by 

. - - 1 nilpotent Lax matrices, which originate from extremal regular (small) D — 4 black holes 

can be obtained as suitable limits of solutions obtained in the diagonalizable case, as we 
show on the generating geodesic (i.e. most general geodesic modulo global symmetries 
of the D = 3 model) corresponding to regular (and small) D = A black holes. As a 
byproduct of our analysis we give the explicit form of the "Wick rotation" connecting the 
orbits of BPS and non-BPS solutions in maximally supersymmetric supergravity and its 
STU truncation. 
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1 Introduction 

The construction and study of black hole solutions in supergravity has a long history. Most 
of the research on this has focused on extremal black holes, not necessarily preserving super- 
symmetry (see e.g. [H El [U H] for reviews). The preservation of supersymmetry makes an 
understanding of the string theory microstates easier, while the vanishing of certain super- 
symmetry variations implies that the second-order equations of motion can be integrated to 
first-order equations, simplifying the explicit construction of solutions. Recently it has been 
shown [Sj El [71 El E] that similar integrations to first-order equations, mimicking the super- 
symmetry variations, can be carried out for some extremal non-super symmetric \10\ [TT\ [T2] 
and even some non-extremal solutions [13\ [H] . 

In this note, we will follow a different approach to solve for black hole solutions in su- 
pergravity, which does not use (hidden) supersymmetry. We use the fact that the black 
hole solutions, after performing a dimensional reduction of the supergravity theory over the 
time direction, are described by geodesic curves on a non-linear sigma model [T3], see also 
|32^ [33l [Ml [35] for early original work. This idea has recently been fleshed out in more detail 
to understand the general structure of BPS solutions [161 [T7], non-BPS attractors [18j, or 
the general properties of extremal and non-extremal solutions [19] . In case the sigma model 
is described by a symmetric space G/H, the geodesic curves are classified in terms of the 
Noether charges of the solution. In particular, when the coset representative L is squared to 
the symmetric coset matrix, Ai = UJ^ 0, then the solution can be compactly written as [15] 

M{t)=MiO)e^\ (1) 

where Q is a matrix containing the Noether charges. This can be seen as a proof of principle 

^L^ denotes the generalised transpose and not the ordinary transpose, see for instance [19j . 
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that the geodesic equations are integrablcl. From the knowledge of Q one can extract already 
a lot of useful information about the black hole solutions [161 fTTl [THl [TU] , but in many cases 
more is needed to understand the physics of the black hole solutions. For instance, the scalar 
U describing the black hole geometry is contained in M and an explicit understanding of the 
black hole geometry requires the full radial dependence of this scalar. However, extracting 
the expressions for the scalars out of the expression for is a hard problem; hence there is 
the need for an explicit integration method on the level of the scalars. 

The integration procedure we will use has first been introduced in supergravity in con- 
structing time-dependent solutions |20 t 1211 [22 | 1231 124j . for which integration through super- 
symmetry is absent from the very beginning. If the solutions possess Killing directions, they 
are described by geodesies on a sigma model of a lower-dimensional supergravity, similar to 
black hole^. The difference between time-dependent solutions and stationary solutions is 
that in the former case the geodesies live in a Riemannian coset, while in the latter case they 
live in a coset of indefinite signature. This difference matters when constructing solutions 
and a main characteristic is that the pseudo-Riemannian case is richer and more involved. 
In the Riemannian case, it was first pointed out in [21] that the geodesic equations could be 
rewritten as a Lax pair and that an explicit integration procedure had been developed by 
mathematicians [25] , This was recently used to construct very non-trivial time-dependent 
solutions in supergravity p3] . 

One of the main purposes of this paper is to demonstrate that this holds also in the 
pseudo-Riemannian case. Namely, the geodesic equations can again be written as a Lax 
pair and the Lax pair is again of a very specific form for which an explicit integration al- 
gorithm has been worked out in the mathematical literature [26] assuming that the initial 
conditions, summarized by giving the Lax matrix at some initial time, are diagonalizable. 
The latter restriction might seem rather unimportant since non-diagonalizable matrices are 
a subset of measure zero in the space of matrices. However, ironically, all BPS states (and 
non-BPS attractive black holes) are described by such initial conditions. It is rather curi- 
ous that exactly non-extremal solutions are easier to describe than extremal solutions in this 
Lax pair approach. This demonstrates that the geodesic approach is orthogonal to using 
(fake) supersymmetry 0. In this paper we will not yet develop the general Lax integration 
algorithm for non-diagonalizable initial conditions; we will however already mention that the 
Lax integration algorithm can be extended to the non-diagonal case. We leave a discussion 
on this issue for an upcoming paper that will contain more technical details |27|. Here, we 
simply want to put forward the principle, and illustrate this with simple examples based on 
the cosets SL(2,M)/ S0(1, 1) and SL(3,M)/ S0(2, 1). A second part of our analysis concerns 
the application to the description of D = 4 black holes. As we shall show, a class of solu- 
tions with non-diagonalizable initial conditions, which can be obtained as limits of solutions 
with diagonalizable Lax matrices, and which therefore can be derived by our algorithm, are 
precisely those which are relevant to the description of non-singular (i.e. having no naked 
singularities) D = 4 black hole solutions. We shall apply our algorithm to the construction of 
the generating geodesic of regular (and small) D = 4 black holes. By generating geodesic we 
mean the solution to the D = 3 model which depends on the least number of parameters such 

■^Proving Liouville integrability of the Hamiltonian system associated with the D = 3 model is a subtler 
issue. It requires the knowledge of a number of conserved quantities in involution equal to the number of scalar 
fields. This problem will be dealt with elsewhere. 

^See [19] for a general unifying explanation of this principle. 

*We refer to [14] for a comparison of the two approaches. 
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that, by applying on it the global symmetries of the model, the most general geodesic can be 
constructed. Such geodesic unfolds in a simpler submanifold of the scalar manifold contain- 
ing a characteristic number of dS'2 factors times 0(1, 1) factors. This solution was originally 
studied in [TH]. Here, we shall further elaborate on it and, as a byproduct, we determine the 
explicit form of the "Wick rotation" connecting the orbits of BPS and non-BPS solutions in 
maximally supersymmetric supergravity and its STU truncation. 

The outline of the paper is as follows. In section [2] we show that the geodesic equations, 
that describe black hole solutions in supergravities with symmetric target spaces, can be 
rewritten in Lax pair form. In section [3] we give a summary of the integration algorithm 
that allows us to integrate the geodesic equations on pseudo-Riemannian symmetric target 
spaces in Lax pair form. This algorithm will be illustrated in section HI using the examples of 
SL(2,M)/S0(1, 1) and SL(3,R)/SO(2, 1). In sectionEl after recalling some results of [19j, we 
shall make general comments on the solutions corresponding to (non-)diagonalizable initial 
conditions and show that the generating geodesic with non-diagonalizable initial conditions 
can be obtained as a limit of solutions with diagonalizable Lax matrices. As an example of 
solutions with non-diagonalizable initial conditions, the generating geodesic corresponding to 
extremal D = A regular (and small) black holes will be explicitly constructed in the maximally 
supersymmetric theory and its STU truncation. Our analysis will show that the Noether 
charge matrices corresponding to such solutions belong, according to their supersymmetry 
properties, to different real sections of the same orbit of the complexification of the global 
symmetry group. Finally, in section [6] we present our conclusions. 

In the final stage of preparation of the present paper, we have learned about the interesting 
paper [28] whose results partially overlap ours. 

2 The geodesic equations in Lax pair form 

As mentioned in the introduction, in describing cosmological and black hole solutions in 
supergravity, one often uses the existence of certain Killing vectors in order to reduce the 
supergravity to a lower dimension. The solutions are then essentially described by geodesies on 
the non-linear sigma model, spanned by the scalar fields of the lower-dimensional theory [15]. 
This non-linear sigma model can be either pseudo-Riemannian or Riemannian, depending on 
whether the reduction to three dimensions includes the time direction or not. In the following, 
we will consider the case in which the non-linear sigma model is a symmetric space G/H , that 
can be either Riemannian or pseudo-Riemannian. In the former case, H is the maximally 
compact subgroup of G, while in the latter case H is a non-compact subgroup of G. In this 
section, we will show that the geodesic equations for the scalar fields can be rewritten in 
Lax pair form, establishing their integrability, irrespective of whether G/H is Riemannian or 
pseudo-Riemannian. The argument proceeds along the same lines as in |21] . where G/H was 
supposed to be Riemannian. 

Let G, EI be the Lie algebras of G and of the (not necessarily compact) isotropy group H 
respectively and let K be the orthogonal complement of H in G, as determined by the Cartan 
decomposition: 



This decomposition is defined through the use of the Cartan involutive automorphism 9, 
which acts as 



G = M + K. 






(3) 
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Since the automorphism preserves the Lie bracket, we have 



[H, M] C M , [M, K] C K , [K, K] C EI . (4) 

Let us illustrate this in case G/H is a maximally non-compact coset, i.e. G is the split real 
form of a complex Lie algebra Gc. In that case, we have 

M = Span{t„} = Span{£;" + ^(E")} , 

IK = SvMKa] = Span{/?,, ^ (S" - d{E''))] , (5) 

where Hi and are the Cartan and the positive step operators (i.e. corresponding to positive 
roots) of the algebra G respectively. For the cosets that originate from a purely space-like 
reduction or from a reduction along the time direction, the action of 6 on the step operators 
is respectively given by [19] : 

purely space-like reduction : 6{E°') = —E~" , 
reduction including time : e{E") = (6) 

where (3o{a) = Hpg{E°') represents the grading of the root a with respect to -f^/jg, the Cartan 
generator that is associated with the internal time directior|f|. Let us now denote by L a 
coset representative. We wish to write the geodesic flow equations in a Lax pair form both 
in the Riemannian and pseudo-Riemannian cases. We shall generically denote by t the affine 
parameter along the geodesic, so that the solution will be described by a suitable dependence 
of the L) = 3 scalar fields on t: cj)^ = 4>^{t) H. The left invariant one-form Q on the coset, 
pulled-back on the geodesic, can be expanded as follows 

n = L-i-^L = <j)^h-^^h = Wta + V^Ka = W + V . (7) 
dt ocp^ 

We identify V as the coset vielbein pulled-back to the one-dimensional space parametrized 
by the affine parameter t. The geodesic action reads 

S = J dtTr{VV) j dtGij{(^)(p^^-^ , (8) 

where the trace is defined in some linear representation of G and where we introduced coor- 
dinates (scalar fields) on the coset via V = KAVf^cj)^ . To compute the equations of motion 
we consider a variation of the action ^ 

5S = 2 j dtTT[V5V] . (9) 



^This generator is normalized so that its adjoint action on G has eigenvalues: (on the generators of the 
four-dimensional isometries and on Hfj^ itself), ±1 (+1 on the shift generators parametrized by the internal 
component of the four-dimensional vector fields and the scalars dual to the _D = 3 vectors; —1 on the shift 
generators associated with the corresponding negative roots); ±2 (on the generators -B±/3o which, together 
with Hpg generate the Ehlers SL(2,R), Ep^ being parametrized by the axion dual to the Kaluza Klein vector 
of the _D = 4 — » _D = 3 reduction). 

^When studying D = 4 spherically symmetric black holes, t will be related to the radial coordinate in the 
Euclidean D = 2> theory originating from a time-like reduction of the four- dimensional one; when studying 
cosmological solutions t will denote the time coordinate in the Lorentzian D = 3 theory arising from D = 4 
through a space-like reduction. 



5 



The variation 5V can be rewritten by using the following identity 

Sn = 6W + 6V = [n, L^^^L] + ^ (L^i^L) . (10) 

The functional L^^^L can then formally be rewritten using the Cartan decomposition 

h-^6h = 6w + 6v (11) 

where 6w and 6v represent the projection of L~^(5L on H, K respectively. Plugging this in 
equation ([TU]) . and projecting the resulting equation onto the IK subspace (using the commu- 
tation relations ^) we find 

6V = [W, 5v] + [V, 5w] + ^5v . (12) 

Finally, upon substituting this expression in the variation of the action and using the cyclicity 
of the trace, we find 

6S = 2 I dtTi 



{[V,W]-^V)6v] , (13) 



from which the Lax pair equation follows 



^V=[V,W]. (14) 

In the following, we will always assume that we work with a coset representative L in solvable 
gauge. As explained in e.g. [201 [ZD [221 [23], this gauge is such that 

W = V>o- y<o , (15) 

where l^>o(<o) denotes the upper-triangular (resp. lower-triangular) part of V Note that 
for pseudo-Riemannian cosets, there is a subtlety in choosing the solvable gauge. This gauge 
can not in general be chosen globally on the manifold for non-compact H. The solvable group 
defined by the Iwasawa decomposition of G with respect to its maximal compact subgroup 
describes local patches of the manifold. Only one of these patches is to be considered as 
physical, namely spanned by the physical fields of the theory. At the boundary of this region 
some of these fields explode, signalling singularities in the corresponding four dimensional 
solution. In the physical solvable patch, time-like and null geodesies (originating from regular 
four dimensional solutions) are complete, while space- like are not, as they reach the boundary 
at a finite value of the " proper time" . We will explicitly show this in section 14.11 where we 
review the simple case in which the scalar manifold G/H \s the two-dimensional de Sitter 
((iS'2) space-time. This example will be particularly instructive since, as will be shown in 
Section [5l the generating geodesic (with respect to the action of G) of regular black holes in 
D = A\s described as a geodesic in a product of dS2 spaces. 

Before starting the discussion of the integration algorithm for the Lax pair equation, let us 
give the expression of the Noether charge matrix in terms of the Lax operator and the coset 
representative. Using ([7]) and ()14p it is straightforward to show that the following matrix 

Q = 2L~^^y'^L^, (16) 



^It is a consequence of Lie's theorem that the generators Ti of the solvable Lie group describing a local patch 
on G/H can be all represented, in a suitable basis, by upper (or lower) triangular matrices. For maximally 
non-compact cosets G/H, the solvability condition can be written as V" = y/^W"'. From ([SJ, © and Q, one 
can then infer that W = V>o — V<o. 
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is a constant of motion. It encodes the conserved charges associated with the invariance with 
respect to the left action of G. Let us stress here that Q is an object of G while V, being in 
IK only transforms under H. The action of a global G-transformation on Q decomposes into 
the action of a global G/H transformation on L whose effect is to move the initial point of 
the geodesic, and the action of a global if-transformation on V. 



3 The Lax algorithm 

In this section, we will consider an algorithm that is useful in solving differential equations 
that can be written in Lax pair form: 

dV 

^ = [V, W] , (17) 

where V and W are x A/^-matrices and W is given in terms of V as 

W = V>o- y<o . (18) 

For the problem of solving the geodesic equations on symmetric spaces, the Lax operator 
V is more specifically defined by 

V{t) = ^ Tr (h-'^KA^ Ka , (19) 

where Ka denote the generators of K. As was shown in section [2l with this definition the Lax 
equation (fTTl) reproduces the geodesic equations on the symmetric space G/H, irrespective 
of whether G/H is Riemannian or pseudo-Riemannian. 

Depending on the symmetry properties of the Lax operator V, algorithms have been 
devised that solve the matrix differential equation (I17p and lead to an explicit t-dependent 
solution for V. After an explicit solution Kioi(^) for the Lax operator V has been found, 
one can generically solve for the scalars that parametrize the coset manifold, by solving the 
following system of differential equations: 

Tr (l-^^Ka^ = TV iV,oi{t)KA) • (20) 

As the left-hand-side of these equations depends on the first derivatives of the scalars, this 
is a system of first-order equations. Depending on the specific parametrization used for the 
coset representative, one can solve this system in an iterative manner. The main difference 
between the Riemannian and the pseudo-Riemannian case lies in the symmetry properties 
of the Lax operator V. For Riemannian cosets, one can choose a matrix representation of 
the Lie algebra, in which all generators Ka are symmetric matrices. The Lax operator V 
is then also given by a symmetric matrix. The algorithm that solves equations (jl7p was 
constructed in [251 129j . For pseudo-Riemannian cosets, the generators Ka are in general no 
longer symmetric. Instead, some of these generators will be symmetric (corresponding to 
the positive signature directions of the coset), while others will be anti-symmetric matrices 
(corresponding to the negative signature directions). Also the Lax operator will therefore no 
longer be a symmetric matrix and will in general be neither symmetric nor anti-symmetric. 
As was shown in [26], the algorithm of [251 [29] can be extended to include also this case. 
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In the following, we will summarize this extended algorithm. The algorithm for Riemannian 
spaces can then be found as a special case of the one outlined below. 

In [26] an integration algorithm for the Lax equations ()17l) (with W given by (jlSp ) is 
outlined, for Lax operators V for which there exists a non-degenerate diagonal matrix S = 
diag(si • • • stv), such that 

V = VS is a symmetric matrix . (21) 

The fact that this integration algorithm allows us to integrate the geodesic equations for 
pseudo-Riemannian cosets is then a result of the following theorem: 

Theorem 3.1. For pseudo-Riemannian symmetric spaces, one can always find a suitable 
space of coset generators lEC, such that there exists a non- degenerate diagonal matrix S with 
the property that 

VKa G IK : KaS = {KASf . (22) 

In order to justify this statement, we note that one can always find a linear representation 
for which the Cartan involution can be written as 

e{X) = -rjX^r, , (23) 

where t] is some diagonal matrix that squares to unity. In general, this matrix r] is given by 

r/ = diag(-lp,lg), (24) 

for some p, q. As 9{]K) = — K, we thus find that 

KaV = {KAvf , (25) 

for Ka S IK, implying that one can take S = rj. Note that also Riemannian spaces obey this 
theorem. Indeed, in that case 5" = H. The formulas that will be given below can then be 
easily adapted to the Riemannian case by taking Si = 1, Vi. 

The algorithm itself is an instance of the inverse scattering method and as such constructs 
the solution Vso\{t) for the Lax operator, starting from the initial conditions contained in the 
Lax operator V at t = 0. The first step in establishing the Lax algorithm consists in realizing 
that the Lax operator can be diagonalized: 

= $A , 

A$ = W^, (26) 
dt ' ^ ^ 

where A = diag(Ai,--- , Xn) is the diagonal matrix containing the eigenvalues of V. The 
matrix $ then contains the eigenvectors of V as its columns. We will denote the eigenvector 
of V with eigenvalue by 



</>(Afc) 



(27) 



\ <pN{h) 
The matrix <I> is thus given by 

CD = [<A(Ai) • • • cPiX^)] = m>^j)]i<i,,<N ■ (28) 
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One can moreover show that ^ can be chosen such that it obeys the foUowing orthogonahty 
relations: 

$5^1$^ = S-^ , <^^S<^ = S . (29) 

Note that in general A can be complex. Also $ will in general be complex, even if A is real. 
Only in the Riemannian case will A and $ be real. 

The t-dependent solutions for the matrix elements of $ are then given by 



^Di{t)D,_i{t) 



SlCll 



</'?(Afc) 



SlCli 



(30) 



The quantities (fP represent the matrix $ at i = 0; i.e. they are obtained from the eigenvalue 
problem 



Fo^o = $0a , 



(31) 



where Vq = V{Q). Note that the eigenvalues contained in A are time-independent, a property 
often denoted as the iso-spectral property of the Lax operator. The quantities Cij{t) in the 
formula ([50]) are given by 



N 



(32) 



k=l 



The Dk{t) are then given by the determinant of the k x /c-matrix with entries SiCij{t), i.e.: 



Dk{t) = det {siCij{t)) 



l<i,j<k 



(33) 



Note that SjCij(O) = 5ij and -Dfc(O) = 1. Furthermore Do^t) = 1. 
The final solution for the Lax operator is then found as: 



N 



[Vsoi{t)]ij = Sj'^s^'^Xk4>i{Xk,t)(l)j{Xk,t) . 



(34) 



fc=i 



4 Examples 

Li this section, we will illustrate the previously outlined algorithm using two examples. The 
choice of the examples is both based on simplicity as well as on physical relevance. The first 
example deals with the SL(2,R)/ S0(1, 1) coset, where the geodesies can be found in a closed 
and rather simple form. We will see that the algorithm indeed leads to the expected results. 
The second example deals with the SL(3,M)/ S0(2, 1) coset. Although still simple, this ex- 
ample is also physically relevant, as it can for instance be used to find black hole solutions in 
four-dimensional Einstein-Maxwell-dilaton theories. This can be done by making use of the 
41? black holes/ 3D SL(3,M)/ S0(2, 1)— geodesies correspondence outlined in the introduc- 
tion. In this paper, we will restrict ourselves to showing how geodesies on SL(3,M)/ S0(2, 1) 
can be produced using the Lax algorithm. The connection between these geodesies and four- 
dimensional black holes will be worked out more explicitly in a forthcoming paper. 
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4.1 The SL(2, M) / S0(1, 1) example 

Let us apply the algorithm to the pseudo-Riemannian coset SL(2, M)/ S0(1, 1). The Cartan 
generator H and positive root E are taken to be 

The Cartan decomposition G = H + IC is then determined by 



Spa.n{l{E + E'^)^ ,K = Span !^H,^ (E - E^) } . (36) 



This is the two-dimensional de Sitter space-time (^6*2) which can be represented by the 
following hyperboloid in M^'^, see Figure [1} 

-{X'^f + {X^f + {Xy = a\ (37) 

where we shall choose = 2. The physical coordinates (solvable coordinates) cp, Xi related 
to the four dimensional fields, span half of this space {physical solvable patch) and are related 
to as follows: e~'^ = X^ + X^ > 0, e~^ x ^ see Figure [TJ These coordinates 
correspond to choosing the coset representative in the solvable subgroup of SL(2,]R) defined 
by the Iwasawa decomposition and will be referred to as the solvable parametrization: 

L(0,x)=e>^^et^. (38) 
The Lax operator in terms of the scalar fields and their time-derivatives read: 

V{t) = - (L"^L-Fr/L^L-^r/) = 

Time-like and light-like geodesies, see for instance |15J, which arise from regular black holes 
in four dimensions, are complete in the physical solvable patch. Space-like geodesies on the 
other hand are not since they cross the boundary X^ + X^ = at a finite value of the affine 
parameter t. These solutions however are related to four dimensional black holes with naked 
singularities. The origin O of the solvable patch, defined by i;^ = x = 0, can be mapped into a 
point O' in the other half of the hyperboloid (X^ + X^ < 0) by the compact transformation: 

W = exp(^ {E - E^)) = J) . (40) 

This transformation, being in the coset, does not leave ry = (-1-1, —1) invariant, but instead 
W~i rjW = —rj. The lower half of the hyperboloid is still described by solvable coordinates 
(p, X, this time parametrizing the coset representative: 

L'((/),x) =e>^^e2^W. (41) 
Let us choose the Lax operator at t = as follows: 
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physical patch 
solvable patches 




\ 



t^to 



Figure 1: The dS2 space-time in M^'^. Some geodesies are represented, for t < 0. The behavior 
for t > is obtained by reflection with respect to the plane containing the axis and passing 
through O. Time-like {\q\ < 1) and light-like (1^1 = 1) geodesies are complete, space-like ones 
are not {\q\ > 1). 



where q = This is a diagonalizable initial condition when \q\ ^ 1. We consider the case 
where \q\ = 1 separately below. The final solution for the Lax operator Vsoi{t) is given by 



where 



where we have defined 



Vu{t) Vi2{t) 

-Vuit) -Vnit) 



Vuit) 



V2 (af{t) - g{t)V^fl^ 



\/2(^cosh VltVa"^ - \/ - k"^ 
— a sinh 



k^ 



fit) 
9it) 



1 _|_ g2V2tVa^-fc2 



1 + e 



(43) 



(44) 



(45) 
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From the above Lax operator, the solutions for the scalar fields can be found. They read, for 
\q\ < 1 as follows 



-0(t) 



X 



X 



Xo + 



-<f>o 



\/l — cosh(x' 

-00 1 + 



2^ 



sinh(x) 

9 

sinh(2 x) 



1 



2(1-9^) 



(cosh(2 x) 



(46) 



(47) 



V2{T^)at, 



where (pQ and xo are the values of (j){t) and x{i) at i = 0. These solutions are regular for 



X < 



log 



, and thus are complete in the limit at 



-oo 



For 



> 1 the 



solutions are obtained from the above expressions by writing x = iy = i \J2 (g^ — 1) at: 



X 



-'t>o 



vV 



Xo + 



-<f>o 



^ — 1 cos(y) — sin(y) , 

1 + J\ sin(2y) 



-4'o 



1 + 



2(9^-1) 



(cos(2y)-l) 



(48) 
(49) 



and are regular for yo < y < yi where yi = Arctan(Y^ — 1) > and y^ = yi — it < ^ 
are the two points in which e~'^ vanishes. One can explicitly check that the above solutions 
dlH]) , (HZl) and p5]) . pU|) do indeed satisfy the geodesic equations. Let us also mention that 
the norm squared of the geodesic is given by 2 (a^ - /c^) . Since the eigenvalues X± of Vq are 



given by Aj 



k ), we see that geodesies with real eigenvalues have positive norm 



squared, while the ones with imaginary eigenvalues correspond to geodesies with negative 
norm squared. 

The behavior of these solutions heavily depends on the values of a and k. For instance, 
upon choosing a = 2 and k = 1, one finds that the eigenvalues of the Lax operator are real. 
The Lax operator is manifestly real and so are the solutions ()46p . ()47p for x and <j): 

p-4>o 



-m 



x{t) 



^/3 



\/3 cosh(\/6t) - 2 sinh(\/6t) 

4^/3 + sinh(2\/6t) 



7-cosh(2V6t) 



(50) 



In the limits t — > ±00, the Lax operator V{t) reduces to a diagonal matrix. One can also 
check that, in flowing from t = —00 to t = +00, the eigenvalues remain constant. In general, 
they will however get permuted on the diagonal during the flow. 

For a = 1 and A; = 2 on the other hand, one finds that the Lax operator has purely 
imaginary eigenvalues A = iy^3/2 and A = — iy^3/2. The Lax operator itself is still real 
however. The solutions for (p and x are given by eqs. (|48p . ()49p . which read: 

1 



x{t) 



xo + 



-00 



=00 



cos{VQt) — cos{VQt) 



V3 



^3 + 2 sin(2^/6t) 
1 + 2 cos(2^/6t) 



(51) 
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The asymptotics for the Lax operator are now very different from the case in which the 
eigenvalues are real. The limits t — > ±00 now no longer lead to a diagonal Lax operator. 
Instead, these limits are not well-defined due to the oscillating character of the functions 
involved. What is physically relevant is the segment of the curve, of finite length, contained 
in the physical solvable patch and parametrized by t G]to, ti[, corresponding to y £]yo, 2/1 [, 
where to < and ti > are the first zeros of e'^ in the neighborhood of the origin O. 

Note that for the Riemannian case, the Lax operator is a symmetric matrix and the 
eigenvalues are always real. The Lax operator thus always reduces to a diagonal matrix at 
t = ±00. Indeed, this phenomenon lies at the heart of the cosmic billiard phenomenon (see 
e.g. [m [231 121]). 

In case the Lax operator at t = is non-diagonalizable, namely \q\ = 1, we can have 
q = e = ±1. The initial value for the Lax operator is: 

where denotes two nilpotent elements such that, if J = ^ (ii^ + E^^) = ^ o) ' have: 

[J, n^] = en^. The solution can be found as the \q\ — > 1^ limit of (j46p . (|47p or, equivalently, 
as the \q\ 1+ limit of (gH]) , (jM]) and reads: 

^-m ^ 6-^0(1- at), (53) 

X{t) = Xo-ee'^ \ (54) 
1 — y Zat 

This solution is regular as long as \plat < 1. The same solution can be found by directly 
solving the Lax equation V + [W, V] = 0. Indeed a nilpotent V{t) in the coset can only have 
the form: V{t) = n^. In this case we will have: W = —e\/2a{t) J. The Lax equation is 

then equivalent to a{t) = — \/2a(t)^, which is solved by a{t) = a/(l — \/2at), where we have 
imposed a(0) = a. Substituting in V{t) and using eq. ([39]) one finds ([53]) and ([M]) . 

This example, analyzed in a certain detail, will be relevant to our discussion of the gener- 
ating geodesic of regular and small D = 4 black holes, which will be done in Section [5] 

4.2 The SL(3, R) / S0(2, 1) example 

4.2.1 Solutions for diagonalizable initial conditions 

Let us now consider the example of geodesies on SL(3, M) / S0(2, 1). We restrict ourselves here 
to diagonalizable initial conditions. The Cartan generators in the fundamental representation 
are given by 

Hi= ^ I , i/2 = I I , (55) 

V 
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whereas the positive roots are given by 



Eu = \ ] , E2S 





£^13 = . (56) 



The subspaces H and IK that determine the Cartan decomposition are then expUcitly given 
by 

M = Span { (^12 + Ef^) , {E23 + E^^) , {E13 - Ef^) } 

K = Span {Hi,H2, {Eu - Ef^) , {E23 - E^^) , (E13 + -^^3) } . 

We will define the coset representative as follows 

Let us illustrate the behaviour of the solutions with diagonalizable initial conditions by giving 
two specific examples. 

All eigenvalues are real: The first initial condition we take, is characterized by the fol- 
lowing Lax operator at t = 0: 



-10 \ 

Vo=\ i(l + 3cosh[7r]) 

^^isMzEl i(l-3coshH) / 



(58) 



The eigenvalues for the Lax operator are all real and given by —1, 2 and —1. Running the 
Lax algorithm, the solutions for the scalars are easily found to be: 

Mt) = Ci + ^V3(4t-log[-l + e6*(-l + cosh[7r]) 
— cosh[7r]]) , 

^log [-1 + e^*(-l + cosh[7r]) - cosh[7r]] , 

^i(V3Ci-C2)^^th [f] 
1 — e6*(— 1 + cosh[7r]) + cosh[7r] 



where Ci, • • • ,6*5 are arbitrary integration constants. The behaviour of the Lax operator is 
very similar to the SL(2, R)-cxample with real eigenvalues. Again, in fiowing from t = —00 
to t = +00 the eigenvalues remain constant and get permuted on the diagonal. This example 
again corresponds to a geodesic with a positive length squared. 





= C2 + 


xo{t) 


= C4 + 


xi{t) 


= C3, 


X2{t) 
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Some eigenvalues are complex: The second initial condition we take, is characterized 
by the Lax operator at t = 0: 

/-I 0\ 
^0= 11- (59) 
V -10/ 

In this case, one eigenvalue is real and given by —1, while the other two eigenvalues are 
complex and given by i (1 + 1^3) and i (l - iVS) E- Upon applying the Lax algorithm, one 
finds the following solutions for the scalar fields: 

Mt) = Ci-^\/3(^-i + log[-3cos[\/3t] + 



VS sin[\/3t] 
lit) = C2 + ^ (st + log [ - 3cos [V^t] + 
VS sin[\/3t]]) , 

2eKv^^i-^2)sin [V^t] 



" ^'^3(-^/3cos[^/3t] +sin[V3t]) ' 
Xi{t) = Cs, 
X2{t) = C5, 

As in the SL(2, M) case with complex eigenvalues, the limits limj^-i-oo V{t) are not well-defined 
and hence do not lead to a diagonal Lax operator, due to the oscillating behaviour of the 
functions involved in the solutions. In this case, the corresponding geodesic is null-like. A 
geodesic with negative norm squared can be obtained by applying the Lax algorithm with 
initial condition 

/-I 0\ 
Vo=\0 12, (60) 
V -2 0/ 

which again has one real eigenvalue and two complex eigenvalues, that are each others complex 
conjugates. 



4.2.2 Regular solutions for nilpotent initial conditions 

After running the Lax algorithm for the initial condition 

/O 1\ 

^0= a , (61) 
V 1 -a J 

*Note that since Vb is real, the complex eigenvalues always come with their complex conjugate. For the 
SL(3, K) example, one will thus always have one real eigenvalue. 
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we obtain the fohowing solutions for the scalars (for 1 — > 0) : 



cosh 



2\/l- aH 



V3 
2 



+ — log[ 



+ cosh[2\/l - aH]] 



2{t) = C2 + Vlog[cosh[2\/l-a2t]] + log[ - 



+ cosh[2\/l -a2t]]j , 

Xo{t) =Ci + a tanh [2^! - aH 

^i(_y3Ci-C2) 



X a+ 



^^^^^ a2 -cosh[2^/^^^] 
^2^i(v^c,+c.)^^ _ ^i(y3C,+C7.)^^ cosh[2yr^t]) , 

X2{t) = C5 + 



-C2 



a e^2/2 _^ e^sinh - aH 



-a? + cosh 



2\/l -aH 



(62) 



In the limit a — > 1, the initial condition ()6ip becomes nilpotent and the solutions (I62|) become 
singular. One can however renormalize the constants Ci and C2 in such a way that the 
renormalized solutions are regular in the limit a ^ 1 and are still valid solutions of the 
geodesic equations in this limit. 

How this renormalization should be performed can be easily seen by looking at the series 
expansion of the solutions around a = 1. For </>i(t) and (p2{t), we get 



Ci + ^log[-l + a] + ^log [-2 - At'] 



+ 



\/3 (3 + 30*2 _^ 40^4^ _ 



+ 0[a-l] 



3/2 



Mt) 



12 + 24*2 

(C2 + ^log[-l + a] + hog[ - 2 - 4*2]) + 
(3 - 18*2 - 56*'^) (a - 1) 



(63) 



12(1 + 2*2) 



+ 0[a - 1] 



3/2 



(Let us ignore for the moment the fact that some terms seem complex instead of real. We 
shall deal with that later.) In both cases, the infinity in taking the limit a ^ 1 comes from 
the second term in the expansion, while all other terms are completely regular. 

The solutions for (/>i(*) and (/>2(*) can thus be made regular in this limit, by renormalizing 
(redefining) Ci and C2 as 



Ci 
C2 



-Vs log 
log - 



Pi 
-1 + a 



P2 

-1 + a 



(64) 
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Note that we will take the new integration constants pi and p2 to be positive. From the series 
expansion above, it is then immediate that the solutions for and 4'2{t) become regular 

in the limit a ^ 1: 

Mt) = ^V3(log[pi(2 + 4t2)]) , 

Mt) = ^(log[p2(2 + 4t2)]) . (65) 

It's an amusing fact that the constants Ci and C2 appear in the solutions for Xoit)i Xi{t) ^^'^ 
X2it) in such a way as to render them regular as well in the limit a ^ 1 after performing the 
renormalization: 

-2V2 {pw2f'H , ^ 

2piP2 + 4plp2i 

Note that the above solutions are real for pi,P2 > 0. One can also check that they obey the 
geodesic equation^. 



5 Relation to the black hole generating geodesic 

As already emphasized, the algorithm discussed in the present note is well defined only for 
solutions corresponding to diagonalizable initial conditions. Although these solutions may 
yield, in certain limits, geodesies with non-diagonalizable initial data, it is not proven that the 
most general solution of the latter type can be obtained in this way. However the generating 
geodesic of regular = 4 black holes (including those with vanishing horizon area, i.e. small 
black holes), defined as the simplest solution capturing all the G-invariant properties of the 
most general one, can be obtained as a singular limit of solutions with diagonalizable Q. 
This generating geodesic was described in [19] as a geodesic within a simple characteristic 
sub manifold oi G/H consisting of a product of dS2 spaces times M+ factors. This analysis 
also applies to extremal solutions, since one can show that the Noether charge matrices of 
light-like geodesies in this simple submanifold have representatives in all the nilpotent orbits 
corresponding to the regular extremal D = 4 solutions (including small black holes). Let 
us first recall the construction made in ^19] of the generating geodesic of regular D = A 
black holes. Then we shall further elaborate on it and work out the explicit transformation 
mapping the BPS and non-BPS extremal solutions. This transformation belongs to the 
complexification G*^ of G but not to G and its action on the generating geodesies is quite 
easy to characterize, having the latter very simple form. For the sake of simplicity we shall 
restrict our analysis to the maximally super symmetric supergravity and to the N = 2 STU 
model. 

We need to first briefly recall the relation between D = A static, single center, asymptoti- 
cally flat black holes and geodesies on G/H. Let us start from a D = 4 extended supergravity 

^We have learned that P. Pre and A. Sorin have independently obtained a similar result. 
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with a symmetric scalar manifold of the form G4,/H4, spanned by scalars (ff and describ- 
ing nv vector fields A^. Being G.i/H^ Riemannian, it is globally described by the solvable 
group iS/i = cxp(Solv4^) defined by the Iwasawa decomposition of G4 with respect to H4 and 
parametrized by (p'": S0IV4 = {(p'^'Sr}- The ansatze for the four dimensional metric and the 
symplectic vector consisting of the electric field strengths and their magnetic duals 
Ga,mi^ oc *SC4/dF^ read: 



ds^ = -e^^dr + e 



-2U 



4 2 



_sinh^ (cr) sinh^ (cr) 

^f,u ' = -e^^C-7W4('A'") -rd^Adr + r sin(g)dgAdyp, (67) 

Ga fivj 2 

where c is the extremality parameter, denotes the 2 ny x 2 ny symplectic invariant ma- 

trix, M.4{(f)'^) = — L4LJ < 0, L4 G S4 being the G4/H4 coset representative in the symplectic 
2 ny-dimensional representation, and = (p^, (/a) the quantized magnetic and electric 
charges. Let us stress that in the present section the radial variable r plays the role of the 
variable t of the previous sections, not to be confused with the time coordinate t. The global 
symmetry group of the D = A theory is G4 which acts simultaneously as an isometry on 
(jf and by means of linear symplectic electric-magnetic transformations on F-^. The charge 
vector r*^ therefore transform in a symplectic representation R of G4. The radial variable r 
is by definition negative. The horizon is located at r ^ —00 while radial infinity corresponds 
to T ^ 0~. 

Upon reduction along the time direction to D = 3 and dualization of the I? = 3 vector 
fields into scalars, we end up with = 3 gravity coupled to a sigma model whose 2 + 2 riy + r),s 
scalars (p^ consist in U and tlic dual a to the Kaluza-Klein vector, coming from the metric, 
the D = A scalars (ff and 2ny scalar fields = (^^, C,h) originating from the vector fields 
and satisfying the relation: = F^. The sigma model metric reads: 

ds^ = GijAcj)^ dcj)-^ = 2{dUf + Grsd(lf dcj)' + ^e-^^ u'^ + e-^^ dZ'^ ■ M4-dZ ,{<6^) 

where uj = da + Z'^ CdZ. The signature of Gjj clearly has 2nv minus signs corresponding to 
the dZ directions. The global symmetry group G of the D = 3 Euclidean theory contains the 
product SL(2,IR)£;xG4, where SL(2, M)^, generated by sl(2,M)£; = Span(ff^p, E|3^^, i?-/3o), see 
footnote 4, is the Ehlers group acting transitively on U, a. Similarly H contains U(1)e x H4, 
where V{1)e is the compact subgroup of SL(2,M)e, generated by Ep^ — E^p^. The charge 
vector transforms in a symplectic representation of U(l)£; x H4 which we shall still denote, 
with an abuse of notation, by R. The solvable group which describes the physical patch of 
G/H is defined by the Iwasawa decomposition of G with respect to its maximal compact 
subgroup (i.e. the compact form of the complexification of H) and its generating solvable Lie 
algebra Solv = Span(s/) is parametrized by (j)^ . The algebra Solv decomposes with respect 
to Solv4, parametrized by (p'^ , as follows: 

Solv = Spsm^Hp^, Ei3g)®Solv4®'R+i, (69) 

where the space R+i = Span(sM) is parametrized by the scalars Z^^ and the grading refers 
to Hpf^. Together with the grading +2 generator E^^, the nilpotent generators sm close a 
Heisenberg algebra: [sm, sn] = ^ 'Cmn Ejj^. The coset representative oi G/H in the physical 
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patch is thus defined as follows: L = exp{aEj3g) exp{^/2 Z^^ sm) exp{(f>^ Sr) exp(C/ H/^g) G 
exp{Solv). 

We can now split the spaces K and M in the Cartan decomposition of G with respect to 
H, as follows: 

K = K2eK4©K(-^^ M = u(1)£;©M4©H(-^\ 

where K2 = Span(//^g, Ef^^ + E_p^-^) is the compact subspace of s[(2, M) £; while IK4 is the non- 
compact space defined by the Cartan decomposition of the algebra G4 generating G4 with 
respect to its maximal compact subalgebra IHI4. The spaces ]K(^) = Span(A;*^) and El(^) = 
Span(sM — sj^^), fc^^ = sm +-5^) transform in the representation R with respect to U(l)£; xi74. 
The former consists of compact matrices and define the negative signature directions on 
G/H. The latter consists of non-compact matrices and generates the Riemannian coset 
H/\[]{1)e X Hii\. We can understand the D = A interpretation of the constants of motion 
encoded in the matrix Q, by restricting, for the sake of simplicity, to geodesies originating at 
radial infinity, from the origin, (f)^ = 0, of the physical patch: L[(/)^(r = 0)] = 1. In this case, 
from eq. ()16p . we can express the corresponding Noether matrix Q = Qo in terms of the Lax 
operator at infinity: Qq = 2 Vg^ G IK, where Vq = V{t = 0). The components of Qo along Hf^g 
and Epg + -E'-/3o are the ADM mass and the NUT charge, the components along K4 are the 
scalar charges while its projection Q^^^ along K^^^ are the electric and magnetic charges: 

q(^) = Qo n ik(^) = -V2r^'CMN k""' = - (/ kA - QA k^) . (70) 

What is the minimal set of generators of K^^) or of H^^^ along which a generic element of 
these spaces can be rotated by means of a \J{1)e x H4 transformation? This minimal set has 
dimension p = rank(///[U(l)£; x//4]) and defined the maximal set of commuting elements { Jfc} 
of E[(^) , or, the maximal number of commuting elements { fcj. } of K^^^ , k = 0, . . . ,p — 1. More 
specifically it defines the dimension of the normal form of the representation R of U(l)£; x H4. 
In other words the minimal number of electric and magnetic charges which characterizes the 
most general geodesic modulo G-transformations (i.e. the generating geodesic) is p. Consider 
for instance the maximal supergravity in D = 4. We have G4 = Ey^y), H4 = SU(8). Upon 
reduction to D = 3 the global symmetry gets enhanced to G = £§(8). In this case H = SO* (16) 
and the D = 3 scalar manifold is therefore Eg(g)/ S0*(16). The quantized charges transform 
in the representation 56 of £7(7), while the central charge matrix Z^b, A,B = 1,...,8, 

and its conjugate Z belong to the 28+i + 28_i of ^{1)e x H4 = U(8). Since we are 
considering geodesies stemming from the origin 0^(0) = 0, at radial infinity the central charges 
{Z^^, Zab)\t=o and the quantized charges are related by a basis transformation. The 
corresponding component of Qo reads: Q^^^ oc (Z^^ kAB — Zab A;^^)|t-=o- It is known that 
a complex 8x8 matrix can always be skew-diagonalized by means of a U(8) transformation. 
If we denote by pk the four real skew-eigenvalues of Z^^ = Zab{t = 0), we can therefore, 
through a suitable U(8) conjugation, bring Q^^^ to its normal form: Q^^^ — > Phk^- In 
this case the four compact generators kk generate the maximal abelian subalgebra of the 56- 
dimensional space K^^^ and p = rank(SO*(16)/U(8)) = 4 coincides indeed with the dimension 
of the normal form. In general we can choose the gener ators (kk) in IfC(^) and (Jfc) in El(^) 
so that, together with Hk = [Jk, kk], generate an SL(2,M)p subgroup of G. In particular 

{Hk, kk) generate a submanifold ^ ^gg'^^^ ^ of G/H where, as we shall see, we can find the 
generating geodesic of extremal (regular) black holes in D = 4. This kind of solutions are 
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generated by a nilpotent Vq (i.e. Q) and moreover can be obtained as a singular limit of 
solutions with diagonalizable Vq. In [TH] it was shown that solutions with diagonalizable Vq 
can be brought, by means of a G transformation, to lie in the following submanifold of G/H: 



/SL(2 



p 



Mn = y-^^^) xSO(l,irP = exp(Kjv), (71) 

where the first p factors are generated by {Hk-, kk). A geodesic in G/H, defined by a diag- 
onalizable Vq and an initial point po = {(Pq) = {(j)^ (0)) ^ G/H, can then be mapped into a 
geodesic on TW^r by first mapping pQ into a point pg of A^at, by means of a transformation in 
G/H, and then using the action of the isotropy group H of p'^ to rotate Vq into the tangent 
space to Mn at Pq. The manifold Mn will be referred to as the normal form oi G/H for 
solutions with diagonalizable Lax operatoJ^. 

Since the p SL(2,M) groups in (I7ip are defined [19j by a maximal set of p-commuting 
generators out of those corresponding to the D = A electric-magnetic charges plus a set of 4 
corresponding Cartan generators, a geodesic oi Mn will describe a four-dimensional dilatonic 
solution coupled to p vector fields. The dilatons parametrizing the S0(1, 1)^~p factor in ([7T]) 
are decoupled from the charges. 

These considerations do not apply to solutions with non-diagonalizable initial conditions, 
such as those originating from extremal, non-rotating, black holes m. D = A, for which 
Tt{V^) = and moreover Vq is nilpotent. In these cases Vq is classified within nilpotent 
orbits with respect to H. In [121 nilpotency condition was obtained for a Vq in the fun- 
damental representation of G, as the Tr(lQ^) limit of a general relation satisfied by the 
non-extremal solutions and reads: 

Vq = (G / £§(8), E8(_24)) ; 

Vi = (G = E8(8) or E8(„24))- (72) 

It turns out that nilpotent elements of the subspace K^v of K have, for different choices of 
their parameters, representatives in all the relevant orbits defined by ()72p . 

In what follows we shall discuss the generating geodesic of regular extremal black holes 
in the maximally super symmetric theory reviewing and extending the analysis in |19j . The 
manifold Mn associated with this model is the same as the one associated with the = 3 
model with scalar manifold: 

G_ ^ SO(4,4) SO(4,4) 

H' SO(2,2) X SO(2,2) S0(4) x S0(4) ' ^ ' 

which originates from the M = 2, D = 4 STU -model coupled to four hypermultiplets. This is 
consistent with the known property that the 5TC/-truncation of the M = 8 model describes 
the generating (seed) solution [9l [30l [3T] of extremal, black holes in the maximal theory. The 
second factor on the right-hand side of (j73p is parametrized by hyperscalars. A generator 
Vq in the [so(4, 4) (so(2,2) ©so(2,2))] space, generating the first factor, transforms in the 
(2,2,2,2) representation of the group SO(2,2) x SO(2,2) = SL(2,M)'' at the denominator. 
It can be written in components as 

^0 = J;Fo"^'^^'^^'^|Ai,A2,A3,A), (74) 



^Note that in [19] the normal form of ^Q^^'^j has been presented in full generality. 



20 



where 1^1,^2)^3)^); = li2, is a basis of matrices; A also labels the supersym- 

metry parameter and is the doublet index of the pseudo-quaternionic structure group 
SL(2,M)o C SL(2,R)^ 0. The matrix Vq in eq. (fTll) can also be seen as an element of a ^ 
q-hit system. Consider regular static black hole solutions. The BPS one corresponds to a 
factorized nilpotent initial condition: y^^i'^2,A3,A _ ^^1,^2,^3 yA g^^^^ yS _ q ^[^^ the 

fundamental of S0(4, 4). In [19] it was shown that the N = 2 non-BPS extremal solution 
with positive quartic invariant ^and vanishing central charge at the horizon) , correspond to a 
different factorization: Vq ^' ^' ^' = C^^'^^'^ v^'' ^ i,j,k = 1,2,3. These solutions, together 
with the N = 2 BPS one correspond to = 8 BPS black hole solutions. Regular solutions 
with non- factorized initial data Vq define non-BPS black holes with negative quartic invariant. 
These correspond to the case in which Vq is an entangled state in the 4 q-bit system. 

Let us now consider the solutions to the smaller model based on the manifold M n ■ Being 
this manifold the product of p = 4 dS2 factors times r — p 0(1, 1) factors, a geodesic on it is 
the product of the p geodesies within each dS'2 factor, which were discussed in Section 14. 1^ 
times r—p geodesies in each 0(1, 1) factor. The dilatons in the S0(1, 1Y~p = S0(1, 1)^ factor 
of A4n, from the N = 2 point of view, are hyperscalars in the N = 2 theory and therefore 
will not be relevant to our discussion of black holes. 

Let us consider extremal geodesies on M.N-, which defined by a nilpotent Vq in K^v n[-So(4, 4)C 
(50(2,2) ©50(2,2))], in the fundamental eight-dimensional representation of S0(4, 4). Being 

nilpotent, Vq can only belong to the generators of the ^ ^30^(2)'' ) factors. 

Recall that the generators {k^, H^) of the f ^57^^") manifold were constructed out of the 



V SO(2) 

normal form of the electric and magnetic charges in D = 4. If the four dimensional STU 
model originates from reduction of a -D = 5 theory, its scalar manifold is described by the 
special coordinate frame and normal forms are the DO — DA — DA — DA (anti-) brane charges 
go, p^ , p^, p^ or the DQ — D2 — D2 — D2 (anti-) brane charges p^ , qi, q2, ^3. Each subspace 
[s[(2,]R) so(l, l)]k within K^r is generated by the nilpotent shift matrices = Hi^ it /c^ = 

1 rx ,.±1 _ 7. _ 1 1^0 1^[T5| 



±1 1) ^^'"^ ih.aX [Jfc,n^] = '^n^, where Jk = 2 \^ oj ^ nilpotent Vq in K^v will 
have the following general form: 

3 

Vo e ©tiK2,K) 050(1, l)]fe , ^0^ = ^ yo = ^n-^" + E^^'^ (75) 

where £k = ±1, k = 0, ... ,3, and the coefficients are related to the four charges of the 
normal form as we will show. The coset representative of I c\f2\ I is the product of p = 4 



SO(2) 

^0 1 
^0 

on the corresponding 2-dim. space. The geodesic with 1^(0) = Vq is the product of p = 4 
geodesies of the form (f^5|) . ([5^ : 



copies of (|38l) : L^r = 0^=0^^'^'''*'°^"^''^'°' where Sk is represented by the matrix 



/to rii 



^'^The supersymmetry variation of the eight fermionic fields in the quarter-maximal theory A'*^''*^'^^ reads: 
^^Each 2x2 matrix, for different values of fe, should be thought of as acting on a different space. 
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where Tik = (1 — V^ik t) and we have chosen 0^(0) = Xfc(O) = 0. Being r < these solutions 
will be regular only for > 0, as we shall assume to be the case. To uplift these solutions 
let us first identify the parameters with D = 4 charges. To this end we write the Noether 
matrix Q = 2 V^: 



Q = V2 ^n^o + ^a^n-^O =^2 ao(i/o + eoA:o) + ^a^(^i-eoA:i) • (77) 




The coefficients V^a^ of are scalar (dilatonic) charges, while the coefficients of kj^ are to 
be identified with the quantized charges according to eq. (f70]l . In particular, we can identify 
go = ^0 clq S'lid = Oj. As for the fields, the relation of the three STU dilatonic fields (pi 
and U to the four dilatons 4>k is: U = \ 0^, = \ (— 0o — 4>i ~\~ 'Pj + (kk), hj: k = 1,2,3 
P^l. The D = A solution then reads: 



This implies that = xo/V2 = and Zi = Xi/V2 = — so that equation Z = 
— e^^C • A^4 • r is satisfied. From eq. ()78p we see that, if all the charges are non- vanishing, 
there is an attractor mechanism at work since at the horizon r — oo, ipi ^^^^^^ , while 

^ = yJAa^ai 02 as , where Ah is the horizon area. The entropy, according to 

the Beckenstein-Hawking formula, reads S = Ah/4: = tt y^Jl^i, where Ia{p,q) is the quartic 
invariant of the representation R of G4 and, on our charges, reads Ii{p,q) = Aqop^ p'^p^ = 
4 £09010203, where e = nfc=o^fc- "^^^ black hole is regular as long as all the four charges 
are non vanishing, in which case one can verify that, by construction, Vq is nilpotent of order 
3 {Vq = 0, Vq / 0) and this is the maximal degree of nilpotency in K]\f. If oq oi 02 03 = 0, 
Vq has a lower degree of nilpotency and the solution describes a small black hole, i.e. a black 
hole with vanishing horizon area. 

We conclude then that the generating geodesic of M = 8 or of the STU model lifts to the 
four parameter dilatonic solution. In what follows we shall focus on the generating geodesic of 
regular solutions, namely for which 13^=0 '^k 0- The fifth parameter of the D = 4 generating 
solution, which is can be identified as with the phase a = Arg(ZZi Z2 Z^^), Z and Zi being 
the STU central and matter charges, can be generated by a G-transformation on the D = 3 
generating geodesic (j78p . Let us now turn to the issue of supersymmetry. The solution is 
supersymmetric if Vq is a SL(2, M)o-doublet. Is this a sufficient condition? From (175]) we 
see that we have 2^ = 16 possible distinct choices for Vq, according to the values of e^- We 
learn from the analysis in [19] that the pseudo-quaternionic structure 51(2, M)o contains the 
generator J = —Jq + Ylf=i -^i- ^0 is an eigen-matrix of J', that is [J', Vq] = ±Vo, only in the 
two cases in which are all equal, for which Vq correspond to the upper or lower component 
of a SL(2, M)o-doublet. These two choices correspond then to BPS solutions (1/2 in = 2, 
1/8 in AA = 8). We have 6 choices for which are not all equal but e = nfc=o^fc ~ 
In these cases, as shown in [19], Vq is no longer eigen-matrix of J but it is eigen-matrix of 
an analogous element on one of the other three SL(2,M) in the isotropy group. The role of 
SL(2, ]R)o is interchanged with one of the remaining three SL(2, R) groups and this corresponds 



^^In the special coordinate frame, where the prepotential has the form — zi Z2 Z3, the three STU complex 
scalars read Zi = bi — ie*^' . 
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to interchange the role of the central charge Z with one of the matter charges Zi . The result 
is a non-BPS solution in the N = 2 model which is still a BPS solution in the M = 8 model 
(since in the latter the four SL(2,M) in the isotropy group are on an equal footing, since they 
form the centralizers of the central charge matrix Zab in the normal form). For these non- 
BPS solutions hip^q) > 0. Finally we have 8 possible choices for which e = nfe=o^'= ~ 
In these cases Vq does not factorize at all and Ii{p,q) < 0. The corresponding solution is 
non-BPS. On the generating geodesic then the only choices which lift to a BPS solution are 
those in which Vq is a SL(2, R)o-doublet. 

We have found then three classes of choices for Vq yielding different kinds of solutions in 
D = 4. These classes are mapped into one another by transformations of the form S^^ Vq S 
whose effect is to switch the grading, n"^* n~^^ , for a subset of the p = 4 terms in eq. 
()75p . keeping the coefficient of n^^ positive. The matrix Sk which does the job on each term 
is Sk = e^^*" giT-'fe since S'^'^ n^'' = n""^*. The matrix S which switches the grading to a 
number i of terms is then: 

e 

s = U^^' C^Q) 

i=l 

If £ is odd, S will cause e = 11^=0 to change sign and thus will map a AA = 8 BPS solution 
Ia{p-,q) > into a non-BPS one Ii{p,q) < 0. Since 5 G G^/G = Ef/Es(s), the orbits of the 
two kind of geodesies will be different real sections of a same E'g'-orbit, consistently with [36] 
. This conclusion is a direct consequence of the analysis in [19], which we have reviewed in 
the present section, about the relation between the generating geodesies of BPS and non-BPS 
black holes. 

6 Conclusions 

In this note we have pursued a different avenue to study the exact integration of all the 
spherically symmetric black hole solutions to supergravity theories with symmetric target 
spaces. This approach applies to all supergravities with more than 8 supercharges and to an 
interesting subset of theories with 8 and less supercharges. Our treatment is not referring 
to (hidden) supersymmetries of the theory, and uses the equivalence between the equations 
governing the radial evolution of the fields in four dimensions, and the geodesic motion of a 
particle on an appropriate pseudo-Riemannian symmetric space. 

We established that the purely t-dependent backgrounds, which reduce to those of a one- 
dimensional sigma model, admit a Lax pair representation and are fully integrable. The 
integration algorithm we have exploited depends more fundamentally on the diagonalizability 
of the initial data. We leave the details of the analysis for extending the algorithm to non- 
diagonalizable initial conditions for an upcoming paper [27] . 

We were able to show explicit analytic formulae for the general integral of simple examples 
like the SL(2, M)/ S0(1, 1) and the SL(3, M)/ S0(2, 1) models. The main message of this note 
is that the integration algorithm is fully explicit and hopefully we have been clear enough to 
illustrate this fact. We have applied this analysis to construct the generating geodesies corre- 
sponding to regular and small D = 4 black holes, which is shown to belong to a submanifold 
consisting of a direct product of dS2 spaces. We have focused on maximal supergravity and 
its STU truncation and, as a byproduct, we have written the explicit form ([7^ of the "Wick 
rotation" mapping the initial data Vq of BPS and non-BPS regular solutions in these models. 
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As a final comment we emphasize that formulating the problem in terms of a Lax pair 
equation is important for ultimately proving the complete Liouville integrability (i.e. the 
global existence of a number of constants of motion in involution equal to the number of 
scalar fields) of the D = 3 model, at least for symmetric spaces. This however is still an open 
issue which we are also working on in |27j . 
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